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Abstract 



o 

j^rj Pushing data traffic from cellular to WiFi is an example of inter-radio access technology (RAT) offloading. 

While this clearly alleviates congestion on the over-loaded cellular network, the ultimate potential of such 



offloading and its effect on overall system performance is not well understood. To address this, we develop a 
general but tractable model that consists of M different RATs, each deploying up to K different classes of 
access points (APs), where each class differs in transmit power, path loss exponent, deployment density and 



\ effective bandwidth. Each class of APs is modeled as an independent Poisson Point Process (PPP), with mobile 

o 

user locations modeled as a further independent PPP, all channels further consisting of i.i.d. Rayleigh fading. 
The distribution of rate over the entire network is then derived for a weighted association strategy, where such 

> 

j~ — ' weights can be tuned to optimize a particular objective. We show that the optimum fraction of traffic offloaded 



to maximize SINR coverage is not in general the same as the one that maximizes rate coverage, defined as the 
fraction of users achieving a given rate. 

I. Introduction 

Wireless networks are facing explosive data demands driven largely by video [1J. While operators 
continue to rely on their (macro) cellular networks to provide wide-area coverage, they are eager to 
find complementary alternatives to ease the pressure, especially in areas where subscriber density is 
high. Complementing the fast evolving heterogeneous cellular networks (HCNs) [|2) with the already 
widely deployed WiFi APs is very attractive to operators and a key aspect of their strategy [3 |. In fact, 
WiFi access points along with femtocells are projected to carry over 60% of all the global data traffic 
by 2015 BH. A future wireless heterogeneous network (HetNet) can be envisioned to have operator- 
deployed macro base stations (BSs) providing a coverage blanket, along with pico BSs, low powered 
user-deployed femtocells [5] and user/operator-deployed low powered WiFi APs. 
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A. Motivation and Related Work 

Aggressively offloading mobile users to WiFi hotspots, however, can lead to degradation of user- 
specific as well as network wide performance. For example, a WiFi network with excellent signal strength 
may suffer from congested backhaul or a certain user's applications may be blocked by a firewall, thus 
reducing the effective rate it can serve at [6]|. Clearly, in such cases any offloading strategy agnostic to 
these conditions is undesirable, which emphasizes the importance of more adaptive offloading strategies. 

Access point association in a multi-RAT setting has been studied in earlier works from both centralized 
as well as decentralized aspect. Fully centralized schemes such as in 0-flU try to maximize a network 
wide utility as a solution to the association optimization problem. Decentralized schemes have been 
studied from game theoretic approaches in lfT0ll - lfl2ll and as heuristic randomized algorithms in |[T3l . 
However most of these works focused on flow level assignment and hence lacked explicit spatial location 
modeling of the APs and users and the corresponding impact on association. The presented work is 
more similar to "cell breathing" lfl4ll - lfT6ll . wherein the BS association regions are expanded or shrunk 
depending on the load. Contemporary cellular standards like LTE use a cell breathing approach to 
address the problem of load balancing in HCNs through cell range expansion (CRE) [2] where users 
are offloaded to smaller cells through an association bias. In particular, this work employs CRE to 
tune the aggressiveness of offloading from one RAT to another in HetNets. Tools from Poisson point 
process theory ifTTl and stochastic geometry |[T8l allow us to quantify the optimal association bias of 
each constituent tier of each RAT which maximizes the fraction of time a typical user in the network 
is served with a rate greater than its minimum rate requirement. 

The metric of rate coverage used in this paper, which signifies the fraction of user population able 
to meet their rate thresholds, captures the inelasticity of traffic such as video services |fT9l whereas 
traditional utility based metrics are more suitable for elastic traffic with no hard rate thresholds. There 
has been considerable effort in the modeling and analysis of HCNs ll2"0ll -|]2"2 l| whereby the location of 
access points of each tier is assumed to be a homogeneous Poisson point process. Furthermore, the 
PPP model for macro cellular deployments has been validated in E3l . [[24]] by comparing with actual 
urban deployments. We take a similar approach where each tier of each RAT is assumed to be a PPP, 
however the key difference from modeling aspect is that the APs of a RAT act as interferers to only the 
user associated with that RAT. Perhaps the most closest work, from modeling perspective, to the work 
presented here is E51 where WiFi hotspots were modeled as a PPP. But the utility function in E51 did 
not consider the effect of the user and AP locations with respect to the association. The present paper 
also characterizes the load - the number of users associated with an AP - which is crucial for finding 
the rate distribution but has not been dealt with in prior work. 
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B. Contributions 

The contributions of this paper can be categorized under two main headings. 

1) Modeling and Analysis. A general M-RAT if-tier HetNet model is proposed with APs of each 
tier of each RAT drawn from a homogeneous PPR This model is validated by comparing the 
analytical results with those of a realistic multi-RAT deployment. 

Association Regions in HetNet: Based on the weighted path loss based association used in this 
work, the tessellation formed by association regions of APs (region served by the AP) is character- 
ized as a general form of the multiplicatively weighted Poisson Voronoi. Much progress has been 
made in modeling association areas of Poisson Voronoi (PV), see G6ll - ll30ll and references therein, 
however that of multiplicatively weighted PV is an open problem. We propose an approximation 
for characterizing the association areas (and hence the load) of an AP which is shown to be quite 
accurate in the context of rate coverage. 

Rate Distribution in HetNet: We derive the rate distribution (CCDF) of a typical user in the 
presented HetNet setting. Rate distribution incorporates congestion in addition to the proximity 
effects that may not be accurately captured by the SINR distribution alone. Under certain plausible 
scenarios the derived expression is in closed form, tractable and thus provides insight into system 
design. 

2) System Design Insights. This work allows the inter-RAT offloading to be seen through the prism 
of association bias wherein the bias can be tuned to suit a network wide objective. 

SINR Coverage: The probability that a randomly located user has SINR greater than an arbitrary 
threshold is called SINR coverage; equivalently this is the CCDF of SINR. In a simplified two 
RAT scenario, e.g cellular and WiFi, it is shown that the optimal amount of traffic to be offloaded 
depends solely on their respective SINR thresholds. The corresponding association bias for max- 
imizing SINR coverage needs to be adaptive to the transmit power of the respective RATs and 
respective SINR thresholds. The maximum SINR coverage under the optimal association bias is 
then shown to be independent of the density of the APs in the network. 

Rate Coverage: We show that in a simplified two RAT scenario the optimal amount of traffic to be 
offloaded for maximizing rate coverage depends on the ratio of the respective resources/bandwidth 
at each RAT and the user's respective rate (QoS) requirements. Specifically, the higher the 
corresponding ratio, the more traffic should be routed through the corresponding RAT. Also, 
unlike SINR, coverage the maximum rate coverage increases monotonically with the density of 
the infrastructure. 



4 



The rest of the paper is organized as follows. Section [TT] describes the wireless HetNet system model. 
The rate coverage expression is derived in Section [III] which is the main technical section of the paper. 
The analysis developed in Section [III] and the PPP model is validated in Section IIII-EI Section ITV] 
explores the system design for optimal offload in a two RAT setting. Numerical results are presented 
in Section |V] with Section [VI] concluding the paper and suggesting future work. 

II. System Model 

The system model in this paper considers up to a K-tier deployment of the access points (APs) for 
each of the M RATs. The set of APs belonging to the same RAT operate in the same spectrum and 
hence do not interfere with the APs of other RATs. The location of the APs of the k th tier of the 
m th RAT are modeled as a 2-D homogeneous Poisson point process (PPP), $ m fc, of density X mk . Also, 
for every pair (m, k) there might be BSs allowing no access (closed access) and thus acting only as 
interferers. For example, subscribers of a particular operator are not able to connect to another operator's 
WiFi APs but receive interference from them. Such APs are modeled as an independent tier (k') with 
PPP $ mfc ' of density \ mk / . The set of all such pairs with non-zero densities in the network is denoted by 
V = [Jm=i with Kn denoting the set of all the tiers of RAT-m, i.e., V m = {{k} : X mk + \ mk > ^ 0}. 
Further, the users in the network are assumed to be distributed according to an independent homogeneous 
Poisson point process $ u with density \ u . 

Every AP of (m, k) transmits with the same transmit power P mk and has an effective bandwidth 
O m k for allocation. The notion of effective bandwidth allows us to capture effects such as backhaul 
congestion and/or firewall blocking stated earlier. In ideal conditions O mk is equivalent to the total 
bandwidth available at an AP of (m, k). The downlink desired and interference signals are assumed to 
experience path loss with a path loss exponent a k for the corresponding tier k. The power received at 
a user from an AP of (m, k) at a distance x is P mk h x x~ ak where h x is the channel gain. The random 
channel gains are modeled as Rayleigh fading with average power of unity, h ~ exp(l). However, 
general fading distributions can be considered at the loss of tractability OTI . The noise is assumed 
additive with power W m corresponding to the m th RAT. Readers can refer to Table ID for quick access 
to the notation used in this paper. 

A. User Association 

For the analysis that follows, let Z mk denote the distance of the typical user from the nearest AP 
of (m, k). In this paper, a general association metric is used in which a mobile user is connected to a 
particular RAT-tier pair if 

(i, j) = arg max T mk Z~ k k , (1) 



5 



TABLE I: Notation Summary 



Notation 


Description 


M 


Total number of RATs in the network 


K 


Total number of tiers of a RAT 


(m, k) 


Pair denoting the fc th tier of the m th RAT 




PPP of the APs of (m, k), PPP of the closed access APs of (m, fc) 


A m fc, A mfc ' 


Density of open access APs of (m, k), density of closed access APs of (m, k) 




Association weight for (m, k); normalized (divided by that of the serving AP) association weight for (m, k). 




Transmit power of APs of (m, k), specifically P m i = 53 dBm, P m 2 = 33 dBm, P m 3 = 23 dBm; normalized 

transmit power of APs of (m, k) 


B m fc; B m fc 


Association bias for (m,k); normalized association bias for (m.k). 


Qfe; at 


Path loss exponent of fc th tier; normalized path loss exponent of k th tier 


w m 


Thermal noise power corresponding to m th RAT 




Effective bandwidth at an AP of (m, k) 


Tmfc 


SINR threshold of user when associated with (m,k) 




Rate threshold of user when associated with (m, k) 


Nmk 


Load (number of users) associated with an AP of (m, k) 


Cmk 


Association area of a typical AP of (m, k) 


S m k', 5 


SINR coverage of a user when associated with (m,k) ; overall SINR coverage of a typical user 




Rate coverage of a user when associated with (m, k) ; overall rate coverage of a typical user 



where T mfc is the association weight for (m, k) and ties are broken arbitrarily. These association weights 
can be tuned to suit a certain network- wide objective. As an example, if T lfe ^> T 2 k, then more traffic 
is routed through RAT-1 as compared to RAT-2. The association weight, T m fc is a general approach, 
special cases include: 

• T mfc = 1: the association is to the nearest base station. 

• T mfe = P mk B mk : is the CRE technique wherein the association is based on the maximum biased 
received power, with B m k denoting the association bias corresponding to (m, k). 

• Further, if {B. mfc } = 1, then the association is based on maximum received power. 
For notational brevity we define 

m A ^ rnk -A PfTtfe -A ^mk i ^ A 
J-mfe = Tp — , -Tmfc — -5 — , E>mk ~ — ana a k — , 

ij ij ij 

which respectively characterize the association weight ratio, transmit power ratio, association bias ratio 
and path loss exponent ratio of a non-serving BS to the serving BS. 

The association model described above leads to the formation of association areas in the Euclidean 
plane as described below. 
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Fig. 1: Association regions of a network with V = {(1, 1); (2, 3)}. The APs of (1,1) are shown as 
hollow circles and those of (2, 3) are shown as solid diamonds. Solid lines show the association regions 
with ^ = 20 dB and dotted lines show the expanded association regions of (2, 3) resulting from the 
use of ^ = 10 dB. 



Definition 1. Association region of an access point is the region of the Euclidean plane in which all 
users are served by the corresponding AP. Mathematically, the association region of an AP of tier 
located at x is 

C Xij = L e R 2 : \\y - x\\ < ' 3 \\y - Z mk \\ a ^W (m, k) e v\ (2) 

The readers familiar with the field of spatial tessellations would recognize that the random tessellation 
formed by the collection {C Xi } of association regions is a general case of the circular Dirichlet 
tessellation lf32l . The circular Dirichlet tessellation (also known as multiplicatively weighted Voronoi 
ll33l Chapter 3]) is the special case of the presented model with equal path loss coefficients which was 
recognized to be the case for HCNs in EOl . Fig. \T\ shows the association regions with two classes of 
APs in the network (V = {(1, 1); (2,3)}, say) for two ratios of association weights ^ = 20 dB and 
^ = 10 dB. The path loss exponent is = 3.5 ("=" is henceforth used to assign the same value to 
a parameter for all classes of APs). 
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B. Resource Allocation 

A saturated resource allocation model is assumed in the downlink of all the APs. This assumption 
implies that each AP always has data to transmit to its associated mobile users and hence users can 
be allocated more rate than their rate thresholds. The available resources at an AP are assumed to be 
shared equally among the associated users. Thus, the rate of a user associated with is given by 

i^ = ^log(l + SINR 4J ), (3) 

where Nij denotes the total number of users served by the AP, henceforth referred to as the load, 
and SINR is the Signal-to-Interference plus noise ratio of the corresponding link. The presented rate 
model captures both the congestion effect (through load) and proximity effect (through SINR). Although 
the above mentioned resource allocation strategy is assumed in the paper, the ensuing analysis can be 
extended to a RAT- specific resource allocation methodology. 

III. Rate Coverage 

This section derives the rate coverage and is the main technical section of the paper. The rate coverage 
is defined as 

K = F(R>p), (4) 

and can be thought of equivalently as (i) the probability that a randomly chosen user can achieve a 
target rate p, (ii) the average fraction of users in the network who at any time achieve rate p, or (iii) 
the average fraction of the network area that is receiving rate greater than the rate threshold p. The rate 
coverage is also exactly the CCDF of rate over the entire network, since the CDF gives F(R < p). 

A. Load Characterization 

This section analyzes the load which is crucial to get a handle on the rate distribution. The ensuing 
analysis uses the notion of typicality which is made rigorous using Palm theory |[T8l Chapter 4]. 

Lemma 1. The load at a typical AP of has the probability generating function (PGF) given by 

G Ni . (z) = E [exp (Xudj (z - 1))], (5) 
where C- V j is the association area of a typical AP of 

Proof: We consider the process $^ U {0} obtained by adding an AP of to the origin of the 
coordinate system which is the typical AP under consideration. This is allowed by Slivnyak's theorem 



[fT8l which states that the properties observed by a typical point of the PPP, $ ij5 is same as those 
observed by node at origin in the process $^ U {0}. The RV Ny denotes the load at this typical AP 
and is the number of users from $ u lying in the association region C Qij of the typical cell constructed 
from the process $^ U {0}. Letting CV, denote the random area of this typical association region, the 
PGF of Nij is given by 

G Nij (z) = E [z N «] = E [exp {\ u C l3 (z - 1))] 

where the property used is that conditioned on CV,, Ny is a Poisson RV with mean X u Cij. ■ 
The following lemma gives the association probability of a typical user to a RAT-tier pair. 

Lemma 2. The probability that a typical user is associated with is given by 

roc ( K \ 

Aij = lixXij J zexp f -n G i:j {k)z 2/&k jdz, (6) 

where 

M 

G ij(k) = ^ mk ^mk k ( 7 ) 



m=l 

If otk = a, then the association probability is simplified to 

X 



^ = „ J . (8) 

3 ZtxGiAk) 



Proof: The result can be proved by a minor modification of Lemma 1 of [21 j. The proof is presented 
in Appendix lAl for completeness. ■ 
As observed from the above lemma, the association probability to (i, j) is directly proportional to 
the corresponding density and association weights. The following two remarks provide an alternative 
interpretation of the association probability. 

Remark 1. The probability that a typical user is associated with the i th RAT is given as Ai = J2j£V t 
This probability is also the average fraction of the traffic offloaded, referred henceforth as traffic offload 
fraction, to the i th RAT. 

Remark 2. Aij is the average fraction of the total area covered by the association regions of the APs 

of (z, j). 

Based on Remark [2] we note that the association area of has mean Below we propose a 
linear scaling based approximation for association areas in HetNets which matches the first moment. 
The results based on the area approximation are validated in Section IIII-EL 



'The term typical and random user are interchangeably used in this paper. 
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Area Approximation: The area of a typical AP of the j th tier of the i th RAT can be approximated 
as 

Cij = C Q 0£) , (9) 
where C c (y) is the area of a typical cell of a PV of density y (a scale parameter). 
Remark 3. The approximation is trivially exact for a single tier, single RAT scenario, i.e. for ||V|| = 1. 



(m,fc)eV ^ mfc 



Remark 4. If T mfc = T and a k = a, then the approximation is exact. In this case, Aij = ^ 
and 

af^)=C £ X mk ). (10) 

V ^ j/ \(m,fc)€V / 

With equal association weights and path loss coefficients the HetNet model becomes the superposition 
of independent PPPs which is again a PPP with density equal to the sum of that of the constituents 
ifTTl and hence the resulting tessellation {C x } is a PV. The right hand side of the above equation is 
equivalent to a typical association area of a PV with density J2(mk)GV ^mk- 

Remark 5. Using the distribution proposed in 11301 for C (y) the distribution of CV, is 

fc tl (c) = f4^V' 5 e X p . (11) 



v i.i 



where T(x) = f °° exp(— t)t x 1 dt is the gamma function. 

To characterize the load at the tagged AP (AP serving a typical mobile user) the implicit area biasing 
needs to be considered and the PGF of the other-apart from the typical-users (N ^j) associated with 
the tagged AP need to be characterized. 

Lemma 3. The PGF of the other users associated with the tagged AP of is 

G No Jz)=3.5^ ( 3 .5 + ^(1-*)) 4 ' 5 . (12) 
Furthermore, the moments of N Qti j are given by 

E fe] = E S(r^)E [C^(l)], (13) 

k=i ^ % i ' 

where S(n, k) are Stirling numbers of the second kind\. 

Proof: See Appendix El ■ 
The moments of the typical association region of a PV of unit density can be computed numerically 
and are available in ||26ll. 



2 The notation of Stirling numbers given by S(n, k) should not be confused with that of SINR coverage, <S. 
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B. SINR Distribution 

The SINR of a typical user associated with an AP of located at y is 

P h U~ a i 

SINRyd/) = v * f ■ w , (14) 

where /i y is the channel gain from the tagged AP located at a distance y, I^k denotes the interference from 
the APs of RAT-z in the tier k. The set of APs contributing to interference are from $jfc\o |J $ iA / V k E Vi 
where o denotes the tagged AP from Thus, 

h,k = Pik Yl h x x~ ak +V lk h x ,x'- ak . (15) 

SINR coverage is defined as the probability that the received SINR is greater than a threshold. The 
spatially averaged SINR coverage of a user when associated with is 

Sijinj) 4 E y pP{SINR ii ( J/ ) > Ti,-}], (16) 

and the overall SINR coverage is 

S= S 'i Ai r (17) 

Interestingly, the distance of a typical user to the tagged AP in Y^, is not only influenced by 

$ij but also by $ mfc V(m, fc) G V. 

Lemma 4. 77ze probability distribution function (PDF), fy^iy), of the distance between a typical 
user and the tagged AP of is 

f Yij (y) = ^//<^p (-7rX>#y/ & 4 (i8) 



k=l 



Proof: See Appendix O ■ 
The following theorem gives the complementary cumulative distribution function (CCDF) of SINR 
over the entire network. 

Theorem 1. The SINR coverage of a typical user is 

S= 2 ^ fVxp ( --3L- - vr r<,-) + y 2 ^ ) dy , (19) 

where 

Dij(k, nj) = P% ak |A ifc Z (jij, a k , TikPik J + A ifc /Z(r ij , a fc , 0)| , 
Gy(k) = jr A mfe t^, Z(a,6,c) = a 2 / 6 f° and SNR 



m=l 



oo 

%3 W 

^ a ' 



1 1 

Proof: See Appendix [Ql ■ 
The result in Theorem \T\ is for the most general case and involves a single numerical integration 
along with a lookup table for Z. In fact, Theorem Q] is equivalent to the earlier derived SINR coverage 
expressions in G3l for M — K — 1 (single tier, single RAT) and that in [12111 for M = 1 (single RAT, 
multiple tiers). 

C. Main Result 

The rate coverage of a typical user in a generic HetNet setting is stated in the following theorem. 

Theorem 2. The rate coverage of a randomly located mobile user in the general HetNet setting of 
Section |77| is 

:3.5r/V, _l A K\ / \ A..\ n / \ A..\-( n+A - 5 ) 



(n + 1) = 2 (n+1)Pi ' /0 " - 1. 



Proof: The rate of a typical user associated with an AP of is given by 

R l3 = ^log(l + SINR u ). (21) 
The probability that the rate requirement of a user associated with is met is 

FiRij > Pa) = P log(l + SINRjj) > ^-J = P(SINR ij > 2«^«/°« - 1) (22) 
= E Nij [S lJ (u ij (N ij ))], (23) 

where Uij(Nij) = 2 Pi 3 N ul°u — 1 and JV i<7 - = 1 + N 0>i j, i.e. the load at the tagged AP equals the typical 
user plus the other users. 

Using Lemma [H ((231) is simplified as 

[<% (iVy ))] = J] ¥(N 0>ij = n)<% (n+1)) (24) 

n>0 

= E— ^r(— J 1, 35 + — J ^ (««(»+!»■ OS 
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Using the law of total probability the rate coverage is 

n = A iA R ij > Pij) (26) 

(i,j)€V n>0 \ J \ V / \ V / 

Using Theorem Q] and Lemma |2] in the above equation gives the desired result. ■ 
Thus, the rate distribution expression for the most general setting requires a single numerical integral 
and use of lookup tables for Z and T. Since both the terms P(iVy = n) and Sy (uij(n + 1)) decay 
rapidly for large n, the summation over n in Theorem |2] can be accurately approximated as a finite 
summation to a sufficiently large value, N max . Particularly iV max = 4A U is found to be sufficient in our 
results presented in Section IIII-EI 

D. Mean Load Approximation 

The rate coverage expression can be further simplified (sacrificing accuracy) if the load at each AP 
of is assumed to equal its mean. 

Corollary 1. Rate coverage with the mean load approximation is given by 

71 = 2 * X v 3/exp J - 7T ( J]) {^"(^ + ^-(fc)} 3/ 2/&fc ) } dj/ , (28) 

{i,j)eV Jo I \k=i J J 

where 

N t3 =E[N ZJ ] = l+ 1 - 28 ^ . (29) 

Proof: Using an approximation for (T23T) with E N .. [«% i u ij{^ij))] ~ «% (*%(E Wij})), Lemma |3] 
gives the first moment of load as E [Ay = 1 + E [N 0jij ] = 1 + ^^E [C®(1)]. Further, using the result 
that E [Cq(1)] = 1.28 11261 . the simplified rate coverage expression is obtained. ■ 
The mean load approximation above simplifies the rate coverage expression by eliminating the 
summation over n. The numerical integral can also be eliminated in certain plausible scenarios given 
in the following corollary. 

Corollary 2. In interference limited scenarios with mean load approximation and with same path loss 
exponents ctk = a, the rate coverage is 

n = V k r ^= r- ( 30 ) 



Proof: For the presented scenario, substituting Wj = and at = a in (1281) simplifies to give the 
desired expression. ■ 
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A careful reader would observe that for a particular rate threshold, the rate coverage is a function 
of the association weights. So it is possible, in principle, to find the optimal association weights and 
hence the optimal fraction of traffic to be offloaded to each RAT so as to maximize the rate coverage. 
This aspect is studied in a special case of two RAT network in Section [IV] 

E. Validation 

In this section the emphasis is on validating the area and mean load approximation proposed in the 
previous sections for rate coverage and on validating the PPP as a suitable AP location model. 

1) Analysis: Through the ensuing results and analysis we aim to verify the validity of the area 
approximation (Theorem [2] and Corollary Q] respectively) and the mean load approximation in the context 
of rate coverage. Consider a two RAT single tier case with V = {(1, 1), (2, 3)}, An = 1 BS/km 2 , 
A 23 = 10 BS/km 2 , X u = 50 users/km 2 , On = 23 = 10 MHz and ct\ = a 3 = 3.5. Figure shows 
the rate distribution obtained through simulation and that from Theorem \T\ and Corollary \T\ for two 
values of association bias. In another scenario, shown in Figure [3l a closed access setting is considered 
where V = {(1, 1), (2, 3), (2, 3')}, with An = 1 BS/km 2 , A 23 = A 23 < = 10 BS/km 2 , A„ = 50 users/km 2 , 
On = 23 = 10 MHz, ax = 3.5 and a 3 = 4. 

As it can be observed from both the plots, the analytic distributions obtained from Theorem \T\ and 
Corollary \T\ are in quite good agreement with the simulated one and thus validate the analysis of the 
previous section. 

2) Spatial Location Model: To simulate a realistic two RAT scenario, the cellular BS location data of 
a major metropolitan used in 11231 is overlaid with that of an actual WiFi deployment 11341 . Along with 
the PPP, a square grid based location model in which the APs for both the RATs are located in square 
lattice (with different densities) is also used in the following comparison. Denoting the macro tier of 
cellular BS as (1, 1) and the femto tier of WiFi APs as (2,3), V = {(1, 1); (2,3)} in this setup. The 
superposition is done such that A 23 = 10 x A n . For the first scenario let B 23 = dB, On = 23 = 10 
MHz, ct\ = a 3 = 3.5 and no noise is considered. Fig. @] shows the rate distribution of a typical user 
obtained from the real data along with that of a square grid based model and that from a PPP, i.e. 
Theorem [2l As shown in the plot, Theorem [2] is quite accurate in the context of rate distribution with 
regards to the actual location data. Fig. [5] shows the rate distribution for a scenario with B 23 = 10 dB 
and «i = a 3 = 4. As per the plot, the accuracy of the PPP is better than that of grid based simulation. 

IV. Design for Optimal offload 

In this section we consider the design of optimal offloading under a specific form of the association 
weight as T m/ t = P m fcB m & corresponding to the biased received power based association or CRE. The 
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objective is to maximize the performance metrics with regards to the association biases {B mfc } and 
hence the traffic offloading fraction. 

For the general settings, the optimum association biases {B m k} for SINR and rate coverage can be 
found using the derived expressions of Theorems \T\ and [2] respectively. We consider below a two RAT 
single tier scenario with q th tier of RAT-1 overlaid with r th tier of RAT-2, i.e. V = {(1, q), (2, r)}. 
Optimal association bias and optimal traffic offload fraction is investigated here in the context of both 
the SIR coverage (i.e. neglecting noise) and rate coverage. 



A. Offloading for Optimal SIR Coverage 

Proposition 1. In interference limited scenarios, with {&k = 1}, A 2r = oAi g and B 2r = bBi q , the value 
of association bias b maximizing SIR coverage is given by 

Big ( Z(rig,a, V 
P 2r \aZ(T 2r ,a,l) 

and the optimum traffic offload fraction to RAT 2 is, 



= ^ : , , (3D 



Z(n g ,a,l) 
Z(r 2r ,a, 1) + Z(n„a, 1) 



Proof: In the described setting SIR coverage can be written as 



5 _ hi (33) 

i ^ ev Y:lx{Di j (k,T ij )+G ij (k)y 



and with V = {(1, q), (2, r)}, A 2r = aX lq and B 2r = bB 



Alq A 2r 



5= ^ „ + 



Xi q Z(r lq , a, 1) + X lq + A 2r (P 2r B 2 r) 2/a A 2r Z(r 2r , a, 1) + A 2r + X lq (P lq B lq ) 2 / a 
1 1 



a(P 2r 6) 2 /« 



Z(r lg , a, 1) + 1 + a(P 2r b) 2 /* Z(r 2r , a, 1) + 1 + a( gJ 6)2/a 

6 opt = arg max j (z(r lq , a, 1) + 1 + a(P 2r 6) 2/a ) + ^Z(r 2r , a, 1) + 1 

_P lg / Z(r lg ,a,l) y /2 
P 2r VaZ(r 2r ,a, 1)/ 

Using Lemma |2l the optimal offloaded traffic fraction is obtained as 

A 2 = X ^ =a L + ( F ^ ) 2/a ) 1 = Z(r lq ,a,l) (34) 

2 G 2r (r) j V p 2r&opt/ Z(r 2r ,a, 1) + Z(ri 9 ,a, 1)' 



The corresponding SIR coverage is 

Z(r 2r ,a, 1) + Z(rig,a, 1) 



Z(r 2r , a, 1) + Z(r lq , a, 1) + Z(r 2r , a, l)Z(r lq , a, 1) ' 



(35) 
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The following observations can be deduced from the above Proposition. 

• The optimal bias for SIR coverage is inversely proportional to the density and transmit power of 
the corresponding RAT. This is intuitive as the denser the second RAT and the higher the transmit 
power of the corresponding APs, the less bias is required to offload traffic. 

• The optimal fraction of traffic/user population to be offloaded to either RAT is independent of the 
density and power and is solely dependent on SIR thresholds. The higher the RAT-1 threshold, r lg , 
compared to that of RAT-2 threshold, r 2r , the more percentage of traffic is offloaded to RAT-2 as 
Z is a monotonically increasing function of r. In fact, if r lq = r 2r , it is optimal to offload half of 
the user population from SIR coverage point of view. 

B. Offloading for Optimal Rate Coverage 

For the design of optimal offloading for rate coverage, the mean load approximation (Corollary Q]) is 
used. 

Proposition 2. In interference limited scenarios, with {&k = 1}, A 2r = a^iq and B 2r = bBx q , the 
association bias b maximizing rate coverage is given by 



arg max <j (z(u lq (N lq ), a, 1) + 1 + a(P 2r 6) 2/a ) 1 + I Z(u 2r (N 2r ), a, 1) + 1 



-'opt 



a(P 2r 6) 2 / c 



(36) 

Proof: The optimum association bias can be found by maximizing the expression obtained from 
Corollary [2] using V = {(1, q), (2, r)}, A 2r = a\i q and B 2r = bB lq . ■ 
Unfortunately, a closed form expression for the optimal bias is not possible in this case as the load 
(and hence the threshold) is dependent on the association bias b. However, the optimal association bias 
b opt for the rate coverage can be found out through a simple linear search using the above Proposition. 
From Proposition [2l it is clear that the optimum offloaded traffic fraction depends on the RAT-specific 
rate threshold of the user and the effective resources 0^. These dependencies are explored further 
in the next section. 

The analysis in this section shows that the optimal traffic offload fraction for SIR coverage and rate 
coverage are considerably different. This is because of the additional dependence of rate coverage on the 
base station loads, rate thresholds and effective resources. An interesting insight from this analysis is that 
SINR coverage and rate coverage exhibit considerably different behavior. For a fixed bias, rate coverage 
always increases with density however for fixed density there is always an optimal traffic offload fraction. 
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This insight might be known by practicing wireless system engineers but here a theoretical makes the 
observation rigorous. 

V. Results and Discussion 

In this section we consider a setting of macro tier of RAT-1 overlaid with a low power tier of RAT- 
2, i.e. V = {(1, 1); (2,3)}. This setting is similar to the widespread use of WiFi APs to offload the 
macro cell traffic. In particular, the effect of association bias and traffic offload fraction on SIR and rate 
coverage is investigated. Thermal noise is ignored in the following results. 

A. SIR coverage 

Fig. [6] shows the variation of SIR coverage with the density of RAT-2 femto APs for different values 
of association bias. The path loss exponent used is a k = 3.5 and the respective SIR thresholds are 
Tn = 2 dB and r 2 3 = 6 dB. It is clear that for any fixed value of association bias, S is sub-optimal for 
all values of densities except for the bias value satisfying Proposition [Q Also shown is the optimum 
SIR coverage (Proposition [U, which is invariant of the density of APs. 

Variation of SIR coverage with the association bias is shown in Fig. [7] for different densities of RAT-2 
femto APs. Again, at all values of association bias, S is sub-optimal for all density values except for 
the optimum density, A opt = (^;) §g|gj} - 

B. Rate Coverage 

Fig. [8] shows the variation of 1Z with the density of RAT-2 femto APs for different values of association 
bias. Variation of rate coverage with the association bias is shown in Fig. [9] for different densities of 
RAT-2 femto APs. In these results, the user density is assumed to be 200 per km 2 , the rate threshold 
p mk = 256 Kbps, the effective bandwidth O mk = 10 MHz and the path loss exponent is a = 3.5. As 
expected, with increasing AP density rate coverage increases because of the decrease in load at each 
AP. The optimum association bias for rate coverage is obtained by a linear search as in Proposition [2l 
however the plots show the exact rate coverage (Theorem [2]). Again, for all values of association bias 
1Z is sub-optimal except for the one given in Proposition [2] It can be seen that the choice of association 
bias can heavily influence rate coverage. 

As noted earlier in Proposition [2} optimum choice of association bias and traffic offload fraction 
depend on the rate threshold (p^) and the effective resources (Oy). To highlight this dependence further, 
the following cases are considered: 

. Case 1: On = 15 MHz, 23 = 5 MHz, Pll = 256 Kbps and p 23 = 512 Kbps. 
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. Case 2: O n = 5 MHz, 23 = 15 MHz, p u = 512 Kbps and p 23 = 256 Kbps. 
The optimum traffic offload fraction and the association bias for these cases are shown in Fig. [TOl As 
can be seen from Fig. [10] for each density value the optimal traffic offload fraction is considerably 
different for the two cases. In particular, larger the ratio of the available resources to the rate threshold 
more is the tendency to be offloaded to the corresponding RAT. 

From the results and discussions presented above it is clear that inter-RAT offload design should 
take into account, not just signal strength and load, but also user's QoS requirements and each RAT's 
bandwidth/resource availability. 



In this paper we presented a tractable model to analyze the effect of offload in a generic M-RAT 
K-tier heterogeneous network setting under a flexible association model. To the best of our knowledge, 
the presented work is the first to study rate coverage in the context of inter-RAT offload. Using biased 
received power based association, it is shown that there exists an optimum percentage of the traffic 
that should be offloaded for maximizing the rate coverage which in turn is dependent on user's QoS 
requirements and the resource condition at each available RAT besides from the received signal power 
and load. Incorporating the effect of traffic on the activity of respective APs, which has been ignored 
in this work, could be an interesting future extension. Although the emphasis of this work has been 
on inter-RAT offload, the framework can also be used to provide insights in inter-tier offload within a 
RAT. Also, the area approximation can be improved further by employing a non-linear approximation. 

Appendix A 
Proof of Lemma [2] 

Proof: If Aij is the association probability of a typical user with RAT-tier pair of then 



VI. Conclusion 




/ 



\ 



(37) 



(m,k)eV 
\(m,k)^(i,j) 

since Z m k denotes the distance to nearest AP in $ m fc. Thus 



(a) 
■ij ~ 



n 




(38) 



(m,k)£V 
(m,k)^(i,j) 




(39) 



(m,k)^(i,j) 
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where (a) follows from the independence of $ m fcV(m, k) 6 V. Now 

P (Z mk >z)=F ($ mfe n b(0, z) = 0) = exp (-vrA m ^ 2 ) 



(40) 



where 6(0, z) is the Euclidean ball of radius z centered at origin. The probability distribution function 

fz mk (z) can be written as 



fz mk (z) = — {1 - P(% > z)} = 27rA^exp(-7rA^ 2 ), \/z > 0. 



Using (ggp, (gOJ) and {41 



^4ij = 27rAij / ^ exp 



z>0 



-IT 



^ ^ Ajfjfc ('I )T(fc) ^ ^ ^ 



(m,fc)eV 



exp(—ir\ijZ 2 )dz, 



which gives ©. 



(41) 



(42) 



Appendix B 
Proof of Lemma [3] 

Proof: As a random user is more likely to lie in a larger association region then in a smaller 
association region, the distribution of the association area of the tagged AP, C[,, is proportional to its 
area and can be written as 



f c > (c) oc cfcJc). 



(43) 



Using the normalization property of the distribution function and (fTTI) the biased area distribution is 



fc'.ic) 



cfcJc) 3.5 3 - 5 A« f\ 



3.5 



exp ( — 3.5^-c 
.Ay 



(44) 



E[Cij] r(3.5)A;\A 3 
The location of the other users in the association region of the tagged AP follows the reduced Palm 
distribution of $„ which is the same as the original distribution since $ u is a PPP lfl"8l Sec. 4.4]. Thus, 
using Lemma Q] and (|44|). the PGF of the other users in the tagged AP is 



GiVo.yO) = E ex P UuCy(z - 1 



.A 



3.5 



A 



-4.5 



exp ( — 3.5— p-c J dc 

"4ij / 

(45) 
(46) 



Using the PGF, the probability mass function can be derived as 

Gfti(0) 3.5 3 - 5 r(n + 4.5) /A U A, 



n! 



r(3.5) 



A; 



3.5 



A 

A,;i 



-(n+4.5) 



(47) 
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For the second half of the proof, we use the property that the moments of a Poisson RV, X 
Pois(A) (say), can be written in terms of Stirling numbers of the second kind, S(n,k), as E [X n 
ELo^S(n,fc).Now 

E W= E [ E [ N o,iMj\_ 
n 



E 



k=0 



^A^S(n,fc)E 



fe=i 



r<'k 



Using (|44|) and the area approximation © 

E W 1 ] _ (A ?J /A,)- (fc+1) E R+^l)] 



E 



and thus 



E fe] = E 



fc=i 



A .4 

Ajj 



(A^/Ai)- 1 !^^!; 



(48) 
(49) 



(50) 



Appendix C 
Proof of Lemma [4] 

Proof: If Yij denotes the distance between the typical user and the tagged AP in then the 
distribution of is the distribution of conditioned on the user being associated with (i, j). Therefore 

¥(Yij > y) = F(Zij > y|user is associated with (51) 

W(Zij > y,user is associated with ^ 
P (user is associated with(z,j)) 

Now 

F (Zij > y, user is associated with = 2n\ij / zexp — ix 

G i:j (k)z 2/&k dz. (53) 

Using (|52]) and ([53]) we get 

¥{Y l3 >y ) = ^l [ Z exp ( -7T Gij{k)z 2/&h \ dz, (54) 
^ Jz> y \ kev k J 

which leads to the PDF of Y^ 

fvM = ^ry^p E GiAk)y 2/ A . (55) 
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Appendix D 
Proof of Theorem CD 

Proof: The SINR coverage of a user associated with an AP of is 

<5; 



. ¥(S\NR, l3 (y)>r, l3 )f Yij (y)dy. 

ly>0 



Now 



P(SINRtf(y) > ni) = P 



> n 



F[h y >y a iP ij - 1 T ij \J2hk + W i 



(56) 



(57) 



E 



exp(-y a %-Pr.i^j i)fc + W i 



(a) 



(58) 



where SNRjj = — and (a) follows from the independence of 1^. The conditional expectation above 
can be written as the moment-generating function (MGF) of the interference, or 



E T . k [exp (-^ P -.%)] = M Iik (//" .-,!',;) . 
Expanding the interference term the MGF is given by 



(59) 



xp | -y^TijPi/Pik <j h * x ~ 



Ctk 



^2 h x >x' ak 



(60) 



(a) 



n M/-. (//' ^i'./'- 

xe$ ik \o 



Oik 



x E 



IK 



(61) 



ik ' \ ■'0 



exp -27rA ifc / U - M fe y^P^-^ idi x exp -2vr V / 1 - M, , ( y Q %-P ifc s'- Qfc x'dx 



exp -27rA ife 



tt 1 + (y^PikTij)-^ 



dx - 2-n\ k i 



x 



1 + (//' I , ,7,' V ) 



dx' 



(62) 
(63) 



where (a) follows from the independence of & ik ,& ik ', h x and (b) is obtained using the PGFL ||T8l 
of $ ife and $ ifc / , and (c) follows by using the MGF of an exponential RV with unity mean. In the above 
expression z ik is the lower bound on distance of the closest interferer in (i, k) which can be obtained 
by using © 



T ijy - a * > T ik z~ a * or z ik > (T ife ) 1A V /a * 



(64) 
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Using change of variables with t = (y aj Pj fc Tjj) 2 / ak x 2 , the integrals can be simplified as 

dx = (i^PifcTfc) 3 ^ / —^-j, (65) 



= (y a tP ik ) 2/a »Z (r i5 , a k , T ttP- 1 ) , (66) 

and 

giving the MGF of interference 

M/ !fc (//• ^IV/) = exp (-vr^P 2 fc /afe {a^Z ( % , a fc) f^T*) + A^Z fa, a fe , 0)}) , (68) 

where 

Z(a,6,c) = a 2 / b ' 



(1)2/6 1 +M 6 / 2 ' 

Using (I58l) . (1591) and <|68]> we get 

P(SINRy (y) > r y ) = exp - vr ^ y 2 / 4 *^ (fc, r tJ ) ) , (69) 

where 



SNR l? 



Dij(k, nj) = P 2 il ak { A ifc Z (nj, a k , P ifc lr T ifc J + A ifc >Z (7^, a fc , 0)| . 
Using (|56l) along with Lemma @] gives 

^ = 2 ^Sy>0 V ^ ( SNP" " ^ ^ {Ai(A: ' ^ + ^ ^ (70) 

Using law of total probability to get 

S= ( 71 ) 

(id)ev 

which gives the SINR coverage of a typical user. ■ 
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Fig. 4: Rate distribution comparison for the three spatial models: real, grid and PPP for a two RAT 
setting with B 2 3 = dB and «i = a 3 = 3.5. 
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Fig. 6: Effect of density of RAT-2 femto APs on SINR coverage. 
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Fig. 7: Effect of association bias for RAT-2 femto APs on SINR coverage. 
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Fig. 8: Effect of density of RAT-2 femto APs on rate coverage. 
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Fig. 10: Effect of user's rate requirements and effective resources on the optimum association bias and 
optimum traffic offload fraction. 



